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Abstract 

We consider the class of autonomous systems x = f{x), where x £ R^", / G C^{T{?"') 
whose phase portrait is a Cartesian product of n two-dimensional centres. We also 
consider perturbations of this system, namely x — f{x) + g{t,x), where g £ C^(R x 
R^") and g is asymptotically small, that is g ^ as t +oo uniformly with respect 
to X. The rate of decrease of g is assumed to be t~P where p > 1. We prove under this 
conditions the existence of bounded solutions of the perturbed system and discuss their 
convergence to solutions of the unperturbed system. This convergence depends on p. 
Moreover, we show that the original unperturbed system may be reduced to the form 
r = 0, = A{r), and taking r G R™, G T", where T" denotes the n-dimensional 
torus, we investigate the more general case of systems whose phase portrait is foliated 
on invariant tori. We notice that integrable Hamiltonian systems are of the same 
nature. We give also several examples, showing that the conditions of our theorems 
cannot be improved. 



1 Introduction 

A. In 1^] the author investigated the asymptotic behaviour of solutions of an R^-system 
X = fix) +g{t,x), 

where the unperturbed system x = f{x) has a phase portrait like a linear centre and the 
perturbation g{t, x) as t ^ +00. In this paper we give natural generalizations of this 
problem. 
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B. First of all we consider an R^"-system 
xi = fi{xi) 

(1.1) 

^ Xn — fni.Xn)t 

where G and fk G C^(R^,R^) for all k = 1, . . . ,n. We also assume that the phase 
portrait of each subsystem Xk = fk{xk) is similar to a linear centre, i.e. origin is a unique 
equilibrium point and the other trajectories arc closed curves surrounding the origin. The 
phase portrait of the whole system (1.1) resembles the Cartesian product 

(R^ X S^) X ■ ■ ■ X (R^ X S^) = R!J: X T", (1.2) 

where R^_ = [0, +oo) and T" denotes the n-dimensional torus. 

The first example of such a system is a linear system x = Ax where all eigenvalues 
of matrix A are pure imaginary and have only prime Jordan blocks. This system can be 
reduced by a linear transformation to a system of n independent mathematical pendulums 

y'i + aiyi = 

y'n + (^nVn = 

where all > 0. The next more general example is a system of n independent Duffing 
equations 

y'i + aiyi = hi{yi,y[) 

y'n + o-nyn = K{yn,y'n) 

where every function hk does not destroy the centre defined by the linear part in each 
equation. We often meet the same systems in different applications of mathematics, me- 
chanics, physics and engineering sciences. There are other important examples of similar 
systems, but we do not mention all of them here. The only exception is a Hamiltonian 
integrable system mentioned below. 
We rewrite system (1.1) in brief form 

x = /(x), (1.1) 

where x = (xi, . . . , x„)^, / = (/i, . . . , fn)^ ■ And we consider a perturbation of (1.1) 

x = f{x)+g{t,x), (1.3) 

where g G C°'^(R| x R^",R^) and function g satisfies the main assumption 

(Al) t*'||g'(f, x)|| uniformly with respect to x on compact subsets ofH?^ as t ^ oo. 

An equivalent formulation of (Al) is: 

(A2) there are positive continuous scalar functions a{s) and j3{t) such that a is increasing, 
P is decreasing, limj_»oo P{t) = and 

tP\\g{t,x)\\<a{\\x\mt)- 
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The proof of this equivalence is standard enough and we omit it here. The foUowing 
theorems are proved for system (1.3). 

Theorem 4.2. If p > 1 then system (1.3) has bounded solutions for t > 0. 

Theorem 4.4. If p > 1 then any bounded solution of (1.3) approaches (orbitally) a 
certain n-dimensional invariant torus of the unperturbed system (1.1). 

Theorem 5.2. If p > 2 then any bounded solution of (1.3) approaches a certain 
solution of unperturbed system (1.1). 

It will be shown that for p = 1 and p = 2 these theorems will no longer be valid. 



C. In Section 3 we construct for the system (1.1) action-angle variables ("quasi-polar" 
coordinates) in which (1.1) has the form r = 0, 6 = A{r) where r = (ri,...,rn)^ are 

"quasi-polar" radii and 6 = {9i, . . . , 9n)~^ are "quasi-polar" angles. 

It gives us possibility to consider instead of (1.1) the more general system in the form 

^ = 

where r G R^, G T", A G C^R."^, Wl \ {0}) (i.e. each component Ak{r) > 0). The 
phase portrait of (1.4) resembles the Cartesian product of a more complicated form than 
(1.2), namely 

(BJ^' X T^) X ■ ■ ■ X (R!^<^ X T"'*) = R![* X T" 

where mi -|- • • • -|- nid = m and ni + ■ ■ ■ + rid = n. The perturbations of (1.4) need to be 
taken in the form 

r = P{t,r,9) 

e = A{r) + Q{t,r,e) ^ ^ ^ 

where P G C^l^gCR^ x R™ x T",R"*), Q G C°;^'g (R^ x R™ x T'*,R") (Let us recall that 
G T" means indeed that 9 G R" and functions P and Q are periodic in each component 

of the vector 9.) 

We shall suppose that P and Q satisfy the following main assumption 

(A3) tP\\P{t,r,9)\\,tP\\Q{t,r,9)\\ ^ uniformly with respect to r on compact subsets of 
11^ and for all 9 ^T'^ as t ^ oo. 

An equivalent formulation of (A3) is: 

(A4) there are positive continuous scalar functions Q:(s) and (3{t) such that a is increasing, 
(5 is decreasing, limj_>(x, /3(i) = and 

t^P{t,rM,nQ{t,rM<oc(\\r\\)m- 



We prove the following theorems for (1.5). 

Theorem 4.1. If p > 1 then system (1.5) has solutions bounded with respect to 
r-coordinates for t > 0. 

Theorem 4.3. If p > 1 then any bounded solution of (1.5) approaches (orbitally) a 
certain invariant torus of the unperturbed system (1.4). 
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Theorem 5.1. li p > 2 then any bounded solution of (1.5) converges to a certain 
solution of the unperturbed system (1.4). 

We note at the end of this section that problems of asymptotic integration are classical 
in the theory of differential equations. If we eliminate the ^-equation and the variable 9 
from (1.5) we obtain a well-investigated in the literature problem of the perturbations of 
zero system (see ||5|- 0], monographs [|| and § give the general theory). From this point 
of view Theorems 4.1-4.4 are quite natural and expected statements. On the contrary, 
Theorems 5.1-5.2 represent more delicate results. We consider the construction of the 
"quasi-polar" variables and proof of Theorems 5.1-5.2 as one of the most important parts 
of our paper. 



2 The existence of a positively-definite integral 

The aim of this section is to prove the following statement. 

Theorem 2.1. Let = fk{xk) be a two-dimensional subsystem of (l-l)- Then it has a 
differentiable positively-definite integral Uk{xk) on R'^. 

This is a known result. But because we shall need further some detailed properties of 
Uk we give an explicit proof of this theorem. 



Proof. For convenience we omit the index k and will write simply x and / instead of x^ 
and fk- We belief that no confusion with the system (1.1) can arise. Let x = (^1,2:2)''' 
and / = {hi, h2)~^ ■ The system x = f{x) can be rewritten in the form 



il = /ll(zi, 2:2), Z2 = h2{zi,Z2). (2.1) 

Consider the system "orthogonal" to (2.1) 

Zi = -h2{zi,Z2), Z2 = hi{zi,Z2). (2.2) 



It is evident that the right-hand vector fields of (2.1) and (2.2) are orthogonal at each 
point. This implies that (2.2) has also a unique equilibrium point at the origin and other 
trajectories of (2.2) cross orthogonally all cycles of (2.1) at a unique point. For if any 
trajectory of (2.2) intersects some cycle of (2.1) at more than one point then there is a 
point on this cycle where right-hand vectors of (2.1) and (2.2) will tangent each other. 

Let ^(t) be a certain nonzero solution of (2.2) and / = (a, b) denote its maximal interval 
of existence. Let P denote the trajectory of ^(t). Without loss of generality we can assume 
that F crosses the cycles of (2.1) from inside to outside with increase of t (otherwise one 
needs to consider the system zi = /i2, Z2 = —hi instead of (2.2)). It follows from Poincare- 
Bendixon theory that ^{t) — > as t — > a+ and ||C(OII ^ oo as t — > 6_. Indeed, system 
(2.2) cannot have closed curves (cycles) because then they must intersect the cycles of 
(2.1) at more than one point. Hence the a-limit set of ^{t) is not empty and may consist 
of only origin. By the same reasons the w-limit set of ^{t) must be empty, i.e. ||C(OII ~^ oo 
as t ^ 6_. It is clear that when time t increases from a to 5 the solution ^{t) crosses each 
cycle of (2.1). We have also that a = — oo because ^(t) remains bounded as t — > a+. 
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Let ip{t,x^) denote the solution of (2.1) with initial condition ip{0,x^) = x^. We define 
the mapping % : \ {0} — > F by the formula 

x{x) = {<f{t, x) -.t eKjnr = <^(r, x) n r. 

It is clear that x € because the cycle <^(R, x) intersects T transversally. Let : F — >■ / 
be the mapping which to every x € L associates the time t such that x = S,{t). Since ^(t) 
is differentiable and ^'{t) 7^ (it is not equilibrium point) we have G C^(r, /). Finally 
we define on R^ \ {0} the function Ui{x) by the equality 

Uiix) = iC'ox){x). 

We point out some properties of C/i. 

1) c/iGCHRnw,/), 

2) uim) = t, 

3) Ui{ip{t,x)) = Ui{x), thus Ui is an integral for system (2.1) on R^ \ {0}, 

4) equation Ui{x) = r defines a unique cycle of (2.1) 

(indeed, curve £,{t) intersects all cycles of (2.1) with increasing t at only one point; 
hence Ui monotonically increases along F; therefore each level curve of Ui{x) consists 
of a unique cycle) . 

5) limjj^o U\{x) = —00. (Indeed, x ^ xi^) (^~^ o x)ix) a = —00). 

By setting U2{x) = exp(Ui{x)) wc construct a continuous positively-definite integral 
for (2.1) on whole R^ (we assume naturally U2{0) = 0). It is clear that C/2(R^) = [0, e''). 
Since differentiability of U2 can be broken only at x = let us consider the behaviour of 
the gradient DU2{x) as a: — > 0. There are two possibilities. 

1) limsup||i:'?72(x)|| < 00. 

2) limsup||£>[/2(x)|| = 00. 

In the first case we define U{x) = ?7|(.t). We thus have liuix^o DU{x) 

= 21ima;^o U2{x)DU2{x) = and U{x) is therefore the required integral for (2.1). 

In the second case we will construct the integral U{x) in the form (/it o U2){x) where 
a scalar function n is chosen to smooth out the discontinuity of DU2 at origin. To do 
this we first define the function z/(s) = sup ||Z)C/2(a^)|| for s G (0,6^*). We have 

V G C(0, e''), since the sets {U2)~^{s) are compact and depend continuously on s. The 
following properties are evident: 

1) v{s) > 0. 

2) lim^^+o i^is) = 00. 



3) i^{U2{x)) > \\DU2{x)\\. 
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Finally we define /x(s) = Jq{1+i'{u)) '^du. It is easy to see that n G Ci[0,e^) and /i(0) = 0. 
We set U{x) = fj,{U2{x)) and verify that U{x) is the required integral for (2.1). All we 
need is to check that DU{x) exists at x = 0. But it is fulfilled, because 



We point out some properties of U . 

1) ?7 G C1(R2, [0,6*)), where h* = e'^^ or /i(e''). 

2) U{x) = <^ X = 0, i.e. U is positively definite. 

3) U monotonically increases along T : dU/dV > 0. 

4) each level curve U{x) = r defines a unique cycle of (2.1). 

The theorem is proved. I 

Note 2.1. For simplicity we will further assume that h* = +oo because it is easy to 
construct a C^'^-function 7 : [0,6*) [0, +00) and consider the new integral ^{U{x)) 
instead of U{x). 

3 The action-angle or quasi-polar coordinates 

In this section we construct new coordinates for more convenient notation of systems (1.1) 
and (1.3). This procedure is a standard one for Hamiltonian systems, see for example the 
textbook 0] . But we would notice that our coordinates are not identical to the hamiltonian 
action-angle coordinates. We introduce them in a more complicated and sly way, so that 
the change of variables does not break the asymptotic properties of the perturbations and 
solutions. 

To each R^-component x^ of vector x we put in correspondence two scalar coordinates 
Tk and Ok- We define the first "quasi-polar" coordinate (the "radius") by = Uk{xk)- 
In accordance to Section 2, the equation Uk{xk) = rk for each > defines a unique 
cycle of Xk = fk{xk)- Moreover the function = Uk{xk) is invertible on F, since U^, 
monotonically increases along F. Let Vk{rk) denote this inversion. We have 



We denote the period of any cycle Uk{xk) = rk by Tk{rk)- It is equal to the time for the 
next intersection with F of the solution passing at t = through the point x^ = Vki'f'k) ■ 
Since F is a tranversal, we have G C^. We now prove the following 



lim \\DU{x)\\ = lim y {U2{x))DU2{x)\\ 




1) ^;fcGCi(R+,r) 



2) Vk{rk) = Tr\{xk: Uk{xk) = rfc}, 



3) Uk{vk{rk)) = rfc. 



Proposition 3.1. The value 1/Tk{rk) is well-defined on R_|_. 
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Proof. Obviously, Tk{rk) for 7^ 0. It remains to prove that rfe(O) 7^ 0. Let us use 
the notations from the proof of Theorem 2.1. Changing in (2.1) coordinates z\ and into 
the standard polar coordinates z\ = pcosa, zi = psina, we get for a 

a = p~^{cos a/i2(pcos a, psin a) — sin a/ii (poos a, psin a)) = p~^F{p, a). 

Obviously, F G C^(R_|_ x R) and F(0, a) = 0. Then Hadamard's lemma implies F{p, a) = 
pG{p,a) where G G C(R+ x R). Hence, a = G{p,a). Let {p{t),a{t)) denote the solution 
passing through the initial point {po,0) at t = 0. Let T{po) be the period of it. It is 
clear that a(T(po)) = 27r. On other hand, a{T{pQ)) = ^ ' G{p{t), a{t)) dt. For any 
po G [0, 1] the solution {p{t),a(t)) belongs to a certain compact set, therefore G{p{t),a{t)) 
is uniformly bounded with respect to po- If T{pq) — > as po ^ 0, then 

rT{po) 

2iT= lim a(r(po)) = lim / G(p(t),a(t)) dt = 0. 
T(po)-o r(po)-o7o 

This contradiction proves the proposition. ■ 
Change time in the system = fk{xk) by formula 
t = Tk{rk)T. 

Since fk = along the solutions we have dt = Tk(rk)dT and 

^ = Tkirk)fkixk). (3.1) 

All solutions of (3.1) have the common period 1, thus the centre defined by (3.1) is an 
isochronous one. Let il^k{T,Xk) be the solution of (3.1) satisfying the initial condition 
"^ki^jXk) = Xk- Wc define the second "quasi-polar" coordinate (the "angle") 6k in such 
a way that ipk{—6k-,Xk) G L. I.e. —6k is the time of the first intersection of solution 
i^ki^^Xk) with r. It is clear that Ok depends smoothly on Xk- Thus the general formula 
for the transformation to the new coordinates is 

Xk = i^k{dk,Vk{rk)). (3.2) 

Let us investigate some properties of the "quasi-polar" variables. For any function 
h{x,y) the notations Dih and denote the partial derivatives of h with respect to the 
first and second argument correspondingly. Let JkifkjOk) be the Jacobian matrix of the 
change of variables (3.2). We have 

Jkin, h) = ^^^^(^'Jq^^^^^ = [D2ipk{dk, Vk{rk)) v'kirk), £>iV'fc(6'fc, Vk{rk))] ■ 

Since i/jk{9k + l,Vk{rk)) = i>k{h,Vk{rk)), the matrix D2tpk{0k,Vk{rk)) is 1-periodic with 
respect to dk- On the other hand D2ipk is the derivative of the solution V'fc with respect 
to initial values, i.e. it is a fundamental matrix of a so-called linear system in variation. 
Hence D2'ipk is nonsingular and D2iJk{0,Vk{rk)) = E. It follows from periodicity that the 
values 

\\D2il^k{0k,Vk{rk))\\, \\{D2il)k{0k,Vk{rk)))~^\\, det D2^pk{0k,Vk{rk)) 
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and det(D2V'fe(^fc) ^fc(^fc))) ^ are bounded away from zero and infinity for each fixed 
and V6'fe G R. 

Since '0ifc(^fe) '"ifc(^fe)) is a solution of (3.1), we have 

Di'ipk{Ok,Vk{rk)) = Tk{rk)fk{ipk{Ok,Vk{rk)))- 

This yields Jk{rk, 0) = [£'-v^(rfc), TA;(rfc)/fc(wfc(rfc))] . But derivative v'j^ is a tangent vector 
to r at the point Vkirk), because Vkifk) is some parametrization of curve T. Hence v'j^{rk) 
is to be orthogonal to the vector fk{vk{rk))- Thus matrix Jfe(rfc, 0) is non singular Vr^ > 0. 

It is easy to check the following statement. Let y = F{y) be an arbitrary smooth 
system and Lp{t, y) be the solution of it with initial condition ip{0, y) = y, then F((p(t, y)) = 
D2ip{t,y)F{y). 

This implies Tk{rk)fk{ipk{6k,Vk{rk))) = D2ipki0k,Vk{rk))Tk{rk)fk{vk{rk)) and therefore 

Jk{rk,dk) = [{D2'>Pk)vk, {D2'ipk)Tk fk] = {D2iJk)[v'k, Tk fk] 
= D2ipk{0k,Vk{rk))Jk{rk,0)- 

Each matrix in the last expression is nonsingular and possesses "good" properties. So 
this equality guarantees that the Jacobian matrix is nonsingular for fixed > and 
arbitrary 9k G R. Moreover, det Jk, det J^^ , \\Jk\\ and H-Z^^H are bounded away from zero 
and infinity. Similar properties of boundedness will also be fulfilled if we take G if 
where K is some compact set of R+ and 6/. as before will be an arbitrary one. 

In the new variables system (3.1) will have the following form 

drk ddu 
ar ar 

The original system Xk = fk{xk) is transformed into 

h = 0, 6k = l/Tk{rk) 

(replacing dt = Tk{rk)d.T). 

Finally the full system (1.1) is transformed into 

r = 0, e = A{r), (3.3) 

where r = (ri,...,r„)^, 9 = {9i,...,9nV, A{r) = (Ai(ri), . . . ,^„(r„))^, Ak{rk) = 

yn{rk). 

Changing the variables in the perturbed system (1.3), we obtain 

G) ^ U)) ^ I'^^-^rMMr))) ^ ^,) . (^<;;;^)) ,3.4) 

(we use natural notations J = diag( Ji, . . . , Jn), ip = {ipi, ■ ■ ■ , tpn)'^ , 

v{r) = {vi{ri), . . . ,'Un('"n))^)- The functions P and Q possess the following properties: 

1) they are in C^'^'q , 

2) they are 1-periodic with respect to each components of 9, 
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3) if (7 satisfies the assumptions (Al) (or (A2)) then P and Q satisfy (A3) (or (A4)). 

Indeed, let K C R" be a compact set and r ^ K. There exist M > such that 
\\{J{r,e))-^\\ < M for all r E K and 9 G R". Hence 

<Mif||5(t,V'(e,^^(r)))|| ^0 as t ^ +00. 

(we note that if r € -fC then cycle ip{9,v{r)) does not leave some compact set in R-^" for 
aU e G R"). 

Thus, we have shown that systems (1.1) and (1.3) may be reduced to systems (1.4) and 
(1.5) while preserving assumption (A3) (or (A4)). 

Note 3.1. We want to underline that it is not so difficult to reduce formally system 
(1.3) to (1.5) (see theory of Hamiltonian system in 0]). The problem is to keep the 
asymptotic properties (A1-A2) as well as the boundedness and covergence of solutions 
under the inverse transformation. For this we use the trick with isochronicity. Without it 
the Jacobian matrix may be unbounded and this manner of reasoning fails. 

4 The existence of bounded solutions 

Let us consider the system (1.5). 

Theorem 4.1. Let P satisfy assumption (A3) or (A4) with p > 1. Then (1.5) has solu- 
tions bounded with respect to r- coordinates. 

Proof. We shall find a comparison equation for ||r||. We have 

j^\\rf = ^(r,r) = 2(r,r) = 2{r, P{t,r,9)), 
here (■, •) denotes inner product. This gives for |[r|| ^ 

|||r|| < \\P{t,r,6)\\ a{\\rmt)/t' < (4.1) 

(we can assume that f3{t) < 1 for sufficiently large t, because /3(t) —>■ 0). Let us consider 
the comparison equation 

j^P = »{p)/t'- (4.2) 

It is well-known (see ||l|) that every function r{t) satisfying inequality (4.1) does not 
exceed the solution p(t) of (4.2) with initial condition /o(io) = lk(^o)||- Thus the existence 
of bounded solutions of (4.2) will imply the same for (1.5). 

Let G{p) = f^l/a{s)ds . Since G'{p) > 0, function G{p) is monotone increasing and 
invertible. The solution of (4.2) with initial condition /o(to) = Po has the form 

G{pit)) - G{po) = {p- ml~P - t'-P). (4.3) 



There are two possibilities. 



Asymptotic Integration of Nonlinear Systems 



207 



A. limp^oo G{p) = +00. Then all solutions of (4.2) must be bounded. For if some 
solution p{t) +CX3 as t — > +00, then the left-hand side of (4.3) tends to co whereas the 
right-hand side tends to {p — l)io~^ < We got a contradiction. 

B. Hmp^oo G{p) = G* < 00. Fix po > and take to so large that G* - G{po) > 
{p—l)tQ~^. Then the solution of (4.2) satisfying the initial condition p{to) = po is bounded. 
Otherwise, tending t +00 in (4.3) and assuming that p{t) — > +00 we get the impossible 
equality G* - G{po) = {p - l)tJ-P. 

Thus the equation (4.2) has bounded solutions in both cases. So the theorem is proved. 



Note 4.1. The case A occurs when P satisfies the well-known inequality (see |Q]) 
||P(t,r,^)|| <M||r||)7(t), 

where J'^j{t)dt < 00 and l/h{s) ds = 00. In this case all solutions of (1.5) are 
bounded with respect to r. Conversely, there are unbounded solutions in the case B as 
the next example shows. 

Example 4.1. Consider the system r = 0, 9 = 1 and its perturbation r = r'^/t'^, 6 = 1. 
The perturbed system has unbounded extendable to infinity solution r = t. All solutions 
lying inside the sector < r < t are bounded and those lying outside it are unbounded 
(moreover, they are unextendable to infinity). 

Note 4.2. The condition p > 1 in Theorem 4.1 cannot be generalized to include the case 
p = 1. See Example 4.2 below. 

Now we give reformulations of Theorem 4.1 and Note 2.1 for system (1.3). 

Theorem 4.2. Suppose that perturbation g in (1.3) satisfies (Al) or (A2) with p > 1 . 
Then (1.3) has hounded solutions. 

Proof. As it follows from Section 3 the system (1.3) written in "quasi-polar" coordinates 
satisfies all conditions of Theorem 4.1. ■ 

Note 4.3. All solutions of (1.3) are bounded if \\g{t, x) \\ < h{\\x\\)^{t), where ^{t) dt < 
00 and l/^(s) ds = 00. 

Example 4.2. The present example shows that condition p> 1 cannot be generalized to 
include the case p = 1. We take the standard two-dimensional centre 

Xi = -X2, X2=Xi, (4.4) 

and perturb it in the following way 

cos t . sin t 

Xi = -X2 + -j— 7, X2 = Xi + ——, (t>l). 
tint tint 

It is evident that cost/ (tint) and s\nt/ {tint) satisfy (A2) with p = 1, a(s) = 1, I3(t) = 
l/\nt. But all solutions of the perturbed system have the form 

xi = (Ci -|- In(lnt)) cost — C2 sint, X2 = (Ci + In(lnt)) cost + C2 sint, 

and are unbounded. 
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Now we return to system (1.5) and show that under the assumptions of Theorem 4.1 
every bounded solution approaches a certain invariant torus of (1.4). 

Theorem 4.3. Let {r{t),9{t)) be a bounded solution of (1.5) with respect to r-coordinates 
and 

iP||P(t,r(t),e(t))|| ^ as t^O (iS) 
with p > \. There exists r* such that 
hm r{t) = r*. 

t— >+oo 

Comment. It means that solution {r{t),6{t)) approaches the invariant torus r = r* of 
(1.4). 

Note 4.4. Evidently, if (A3) or (A4) holds then every bounded solution of (1.5) satisfies 
{A3). 

Proof. We have 

/•t j-t 

\ ds. 



it)=r{to)+ [ r{s)ds = r{to)+ [ P{s,r{s),9{s)) , 

Jtn Jtn 



'to J to 

Assumption (^3) implies that \\P{s,r{s),6{s))\\ ds converges. Hence the limit 
lim r{t)=r{to)+ / P{s,r{s),e{s)) ds = r* 

Jto 

exists. The theorem is proved. I 

Note 4.5. The condition p > 1 cannot be generalized to include the case p = 1, see 
Example 4.4. 

Note 4.6. The next example demonstrates that assumptions of Theorem 4.3 don't imply 
even orbital convergence of bounded solutions of (1.5) to certain solution of (1.4). However 
this convergence occurs if € (see [^), see also Note 5.2. 

Example 4.3. Consider the following system 

ri=0, r2 = 0, ei=ri, ^2 = ^2, (4.5) 
and its perturbation 

h = l/t\ r2 = 2/t^ ^i=fi, ^2 = r2, (4.6) 
Solutions of (4.5) are 

ri=r?, r2 = rl = r?t + 0?, 02 = r^^t + Ol (4.7) 
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where r^, 9^, 62 are arbitrary constants (initial values). Solutions of (4.6) are 

fi=f?-l/t, f2 = f^-2/i, ei = fit-lnt + e^, 02 = fit - 2lnt + 9^. (4.8) 

All solutions of (4.6) arc bounded with respect to fi and f2- It is clear that fi and 
f2 f2 as t 00. The orbital convergence of (4.8) to (4.7) means that = k = 1,2 
and there is an infinitely monotonically increasing function s{t) such that 

lim (Oiit) - eMt))) = lim (^2(t) - e2ism = 0, 

t— >oo c— >oo 



or 

lim (r?i - Int + - r?s(t) - e?) = 0, (4.9) 

t—>oo 

lim {r^t -2lnt + 9^- r^s{t) - 9^) = 0. (4.10) 

t— >oo 

Prom (4.9) we have s{t) = t - (Int -9^ + 9^)/r^ + o(l). Substituting this in (4.10) we 
obtain 



hm ((^ - 2) Int + 9l - 9^ + (0"? - ^?)/r? + o(l)) = 0. 

However, this is possible if and only if = 2r^, 9^ = 9^, ^2 = ^2- Consequently there are 
solutions of (4.6) which don't orbitally converge to any solution of (4.5). 

Now we reformulate Theorem 4.3 for system (1.3). 

Theorem 4.4. Let x{t) be a bounded solution of (1-3) and 

tP\\g{t,x{t))\\ ^0 as t^ +00 
with p > 1. Then x{t) approaches a certain invariant torus of unperturbed system (l-l)- 
Proof. It follows immediately from Theorem 4.3 and Section 3. ■ 

Note 4.7. For Theorem 4.4 Notes 4.4-4.6 are also valid. 

Example 4.4. The next example demonstrates that assumption p > 1 cannot be gener- 
alized to include the case p = 1. Let's consider again the two-dimensional linear centre 
(4.4) with the following perturbation 

cos(ln(lnt)) . sin(ln(lnt)) . ^, 

Xl = -X2 + W " xu X2 = Xi+ \} " x2, i>l. 4.11 

tmt tint 

Obviously every bounded solution of (4.11) satisfies the conditions of Theorem 4.4 with 
p=l. But (4.11) has the bounded solutions 

xi = rocosiexp(sin(ln(lni))), X2 = ro sint exp(sin(ln(lni))), 

which oscillate between the circles with radii roe and roe~^. Hence no one of them can 
approach any circle x\ = r cos t, X2 = r sin t, which is just the solution of the unperturbed 
system. 
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5 The convergence of bounded solutions 

We begin with the main statement for system (1.5). 

Theorem 5.1. Let (A3) or (A4) hold with p > 2. Then each bounded solution of (1-5) 
converges to a certain solution of the unperturbed system (1-4) ■ 

Proof. Let {r{t),9{t)) denote any bounded solution of (1.5). Theorem 4.1 imphes the 

convergence r(t) ^ r* as t ^ +oo. Every solution of (1.4) lying on the invariant torus 
r = r* has the form {r*,A(r*)t + ^o)) where is the initial value. It is necessary and 
sufficient to prove the existence of a such that 

lim {9{t) - A{r*)t - Oq) = 0. 
Thus it is sufficient to prove that the limit value 9o = \imt-^+c>o{9{t) — A{r*)t) exists. Since 

9{t) - A{r*)t = 9{1) - A{r*) + J\Hs) - A{r*)) ds 

= 9(1) - A{r*) + j\A{r{s)) + Q{s, r{s),9{s)) - A{r*)) ds, 
it is sufficient to prove that the obtained integral converges. By (A4) we have 

\{s,r{s),9{sm<a{r{s))P{s)/sP. (5.1) 
There is a constant M such that a{r{s))P{s) < M for s > 1. Hence 

/oo /"OO 
Q{s,r{s),9{s))ds\\ < M j l/s^ ds = M/{p - 1). 

Next we need to estimate ||^(r(s)) — ^(r*)||. First of all we notice that 

f{u)du\\< / \\P{u,r{u),9{u)) du\\ 

.-oo Js 
r+oo r+oo 

< a{r{u))(3{u)/uPdu<M j l/u^ du = Ms^'^ /{p - I). (5.2) 

J s J s 

Since A G C^, it satisfies the Lipschitz inequality with some constant L (on such a compact 
set K from R!}? that r(s) G iiT for s > 1 ) 

\\A{r)-A{r*)\\<L\\r-r*\\. (5.3) 

It follows from (5.2) and (5.3) that \\A{r{s)) - A{r*)\\ = 0{s^~p). Inequality j3 > 2 implies 
convergence of \\A{r{s)) — A{r*)\\ ds . Hence limit 

poo 

lim {9{t) - A{r*)t) = 9{1) - A{r*) + / (^(r(s)) + Q{s, r{s), 0(s)) - A{r*)) ds 

t^+oo 

exists. This proves the theorem. ■ 
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Now we give the formulation of the analogous theorem for (1.3). 

Theorem 5.2. Let (Al) or (A2) hold with p > 2. Then each bounded solution of (1-3) 
converges to a certain solution of the unperturbed system (1.1). 

Note 5.1. The condition p > 2 cannot be generalized to include the case p = 2 as the 
following example shows. 

Example 5.1. Consider the system 

xi = -xl, ±2 = xl, (5.4) 
and its perturbation 

3 1 + lnt . 3 1 + lnt , , 

^^ = -^^+(U^^^' ^^ = ^^+(iW^^- ^'-'^ 

The cycles of (5.4) arc xf + = C, C € R+. Let (Cnt, Snt) denote the solution of 
(5.4) passing through point (1,0) as t = 0. Changing the variables by the formulas xi = 
rCnO, X2 = rSnO (action-angle coordinates) we transform (5.4) into 

f = 0, e = r'^, (5.6) 

and (5.5) into 

The solutions of (5.6) are r{t) = tq, 9{t) = + f^t and those of (5.7) are f{t) = 
roexp[-l/(tlnt)], 6{t) = 6*0 + rg /2 exp[-2/(s In s)] ds. Obviously, limf(i) = ro. But 
lim {9{t) — A{ro)t) from proof of Theorem 5.1 does not exist. Indeed, 

t— >+oo 

lim {d{t) - A{ro)t) = lim (6'(i) - rgi) = lim rg( f exp[-2/(s In s)] ds - t) 

t— >+oo t— >+oo t^+co J2 

= lim rli [ (exp[-2/(slns)l - l)o?s - 2) = -oo, 

t-^+oo J2 

because obtained integral diverges. This follows from the asymptotic equivalence 

exp[-2/(slns)] - 1 2/(slns) 

and divergence of 2/{slns) ds. 

Note 5.2 (about convergency in Section 4). We wouls like to make an interesting 
remark about the possibility to decrease the lower bound in the inequality p > 2 in 
Theorems 5.1 and 5.2. Suppose A{r) = is a constant. For system (1.1) it is equivalent 
to Tk{rk) = Tfc. It means that each subsystem Xk = fk{xk) in (l-l); k = 1, . . . ,n, defines 
an isochronous centre. For (1.5) we can also talk about a property of isochronicity. Then 
the statements of Theorems 5.1-5.2 are valid for p> 1. Indeed, we have in this case 

/•oo 

lim {e{t) - A{r*)t) = 9(1) - A{r*) + / Q{s, r{s), 9{s)) ds, 
where convergence of this integral follows from (5.1). 
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The property of isochronicity in Hamiltonian systems (and not only for such systems) 
is discussed in numerous papers. We refer only to where further references can be 
found. Only for asymptotic perturbations of such systems, our Note 5.2 is valid. 

We note that the main reason of appearance of the condition p > 2 is just the necessity 
to estimate A{r{t)) — A{r*). Hence we may suppose that Q satisfies (A3) with p > 1 (but 
P with p > 2 as before) and Theorem 5.1 will also be valid. 

Finally we point out a possibility of the situation for system (1.1) when some subsys- 
tems of it are isochronous centres and all other not. It is not difficult to see that for 
"isochronous" angles 9k the orbital convergence (Note: even usual convergence) is realized 
in Theorem 4.4. 
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